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Abstract 

Recently it has been suggested that junctions between materials 
with different parity violating properties would be characterized by 
diffusion layers, analogous to those in the p-n j unction This 
remark is amplified by a fuller investigation of two related parity vi- 
olating effective Lagrangians, which possess a kind of duality. It is 
shown that gauge invariance and energy conservation are sufficient 
to determine the behaviour at the interface. This leads to modifi- 
cations of normal parity- violating electrodynamics. The coupling of 
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an interface to an external system is a natural solution to the de- 
ficiencies of Maxwell-Chern-Simons theory. A heuristic model of a 
transistor-like device is discussed which relates to recent experiments 
in device technology. Radiative corrections to Chern-Simons theory 
induce a local magnetic moment interaction whose lagrangian is every- 
where gauge invariant. The effects of this interaction are compared to 
Maxwell-Chern-Simons theory. The dispersion of classical waves for 
these models is computed and the laws of reflection and refraction are 
found to hold despite the lack of P and T invariance. The magnetic 
moment dispersion is gapless in contrast to the Chern-Simons disper- 
sion except in the case of a scalar field which is covariantly constant. 
Both models exhibit optical activity (Faraday effect). 

1 Introduction 

The Chern-Simons term is widely exploited in the construction of effective 
theories breaking parity and time-reversal invariance. Although originally 
the term was introduced as a way of providing the Yang-Mills gauge field 
with a gauge invariant mass[|l], H, H, it has since appeared more often in 
the condensed matter literature in discussions of systems like the fractional 
quantum Hall effect and anyon sup er conduct ivityH]. These models postulate 
the Chern-Simons term as an effective action for an unknown microscopic 
theory, the coefficient of which takes on a constant value which is chosen or, 
in principle, determined from the underlying physics. 

In this paper we explore more carefully the possible roles of the Chern- 
Simons Lagrangian by investigating systems in which its coefficient - the 
physical potency and sign of the parity violating effects - takes on different 
values in different parts of the system. An additional model which couples 
a gauge invariant current to the dual of the field strength is also considered 
and compared to the usual Chern-Simons term. This model reduces to the 
Chern-Simons expression in a special case. 

The junction scenario described in this paper, although motivated on 
theoretical grounds, could have some bearing on recent experiments in which 
participating electrons are characterized by predominantly a single spin direction P, 
p. Redlich has shown that spin polarized Dirac fermions give rise to an ef- 
fective Chern-Simons theory when quantum corrections are accounted for[0 
and thus a Chern-Simons theory describes such a two-dimensional electron 
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gas. 

We begin by formulating the simplest junction in terms of an action 
principle. This consists of two regions in a (2+l)-dimensional space, which 
meet on the line Xi — 0. The two regions are considered to have physically 
disparate properties so that a boundary condition is implied for the physical 
fields in addition to the relevant field equations on either side of the boundary. 

2 Formalism 

To illustrate the variational formalism we are using, consider the simplest 
case of plain electromagnetism at a material junction. In Maxwell theory, the 
only variables in a physical medium are the dielectric properties of matter; 
it is sufficient to describe these properties in terms of a conserved electric 
current variable Jtot{^) whose value is position dependent and a dipole 
field P^i, = D^i, — F^i,. The action is given by 

S^J dV,{~D^^F^, - JV^J (1) 

where F^^ = d/^A^ - dyA^ {F^^ = B, = E'), dV^ = dtdxidx2 and the 

metric signature is 77^1^ = (H ). Units are chosen such that h — c — 

^0 = A*o = 1 a-nd latin indices refer to spatial dimensions only. We shall have 
occasion to consider both sharp and soft interfaces. To take account of the 
polarization of media, one must allow the possibility that currents will run 
along the surface of a sharp interface at xi = 0. One can, with a certain 
freedom, either include this effect in P^^ or introduce the current explicitly 
as a surface current. In the following we include both, anticipating the work 
in later sections. Since we are modelling the interface as an abrupt change 
in the physical properties of our system, these currents will be proportional 
to a delta-function in the xi direction. We write 

J^OT - J' + r S{x^) ^ + (2) 

— is expected from the geometry (this is not the case in the Chern- 
Simons theory) and the total current will be conserved dfj^Jj^Qrp — 0. Varying 
the action and integrating by parts leads to 

5S ^ j dV^lSA" d^{F^, + P^y) - J^5A^} 
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+ J dV^d^ [-5A'' {F,, + P^,)] - 5A^j, = (3) 

Requiring that Q leads to 

= J dV^{5A^ d^D^, - r5A^} 

+ j dtdx2 dxi{d^ [-SA" Di^] - dA^j^ 6{xi)} = 0. (4) 

Moreover, if the regular field equations (in the absence of a boundary) are 
to apply arbitrarily close to Xi = 0, then the boundary condition 

SSboundary = 1™ J dtdX2 J^^ dXi{d^ [-6 A" {F^^ + Pi^)] - dA"]^ 6{Xi)} = 

(5) 

is implied. Here it is assumed that the components of the vector potential A'^ 
are continuous across the boundary. This assumption may be unnecessarily 
restrictive if one considers the case in which a contact potential characterizes 
the junction and will be relaxed later. This leads to the form 

A(Fi^ + Pi^)+j^ = fi = 0,2 (6) 

where A means the change in value across the boundary. The field equations 
outside of the boundary are given by 

d^D^, = J, (7) 

Clearly one can absorb the contribution due to surface currents into the more 
usual form of a polarization tensor. The surface polarization tensor may be 
defined by 

d'P^u = Jt (8) 

which is added to P^^. The antisymmetric components of the total polariza- 
tion are Pi^ = —M, P^^ = P^, enabling to be written in the standard 
form 

AD^ = A{E'^ + P^) = 
AH = A{B - M) = 0. (9) 

where P^ and M are the polarization and magnetization respectively. It 
also follows immediately from the assumption about the continuity of the 
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vector potential that /S.E2 = A{d2AQ — 99^2) = 0. The conservation of the 
total current gives us a relation between and the change in across the 
boundary (AJ^). 

d''{J^+j^5ix,)) = (10) 

Integrating this equation with respect to Xi from — e to +e and consider the 
limit e — * 0, one obtains 

A Ji + d^j, = 0, (ji = 0) (11) 

The main point to note from this exercise is that it is unnecessary to 
explicitly introduce fields D and H to account for polarization and magne- 
tization effects specifically on the boundary: it is sufficient to include the 
possibility of surface currents. We shall therefore not refer to these fields 
again. 

3 Maxwell- C hern- Simons Theory 
3.1 Gauge invariance at a junction 

The action formalism does not yield an explicitly gauge invariant boundary 
condition when the Chern-Simons term is considered. The Abelian Chern- 
Simons term is gauge invariant only in the absence of boundaries, being 
quadratic in but only linear in the derivatives. Consider the action, 

S = I dVA-lp^'F^, + ^e^-^A^d.Ax - J^oT^^} (12) 

The variation of this action, without further stipulation, leads to the field 
equations and associated boundary condition 

^^F^, + ^^e,px^PA^ = J, (13) 
A[-FiA + ^/ie^iAA^']-jA = (14) 

Gauge invariance of the boundary condition implies a restriction on A^. One 
is therefore led to consider the effect of the gauge transformation A'^ = 
A^ — d^O on the action (|T2]) over a region — e < Xi < e as e — > 0. 

5S = S[A'^]-S[A,] = J dtdx2 J^^dxi{-d'[^fie,ixd-A^+Ji]-{d^j^)6ixi)}9 

(15) 
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where ji = 0, jo, 72 are independent of xi and = If the gauge choice 
6 is to be arbitrary one has 

A(^/iE2 + Ji) + 5^> = 0. (16) 

On use of ( pi]) this simply becomes 

AfiE2 = (17) 

since the continuity of the vector potential implies that AE2 = 0. This result 
indicates that the only valid boundary condition between Chern-Simons me- 
dia with different coefficients is one in which £"2 = on the boundary. This 
boundary condition does not permit the passage of electromagnetic waves or 
information. The vanishing of E2 implies the form 

= 9,e(x) ^ + 1 (18) 

for some scalar field ^[x). This solution may be used in the action as a 
restriction on the allowed variations. The analogue of equation (^) is then 

?>S^oundary= J dt dX2 J^^ dxi{d' d" F^, - ^d,{fie^ixd^ A^) + d'' J ^6{x,)}6^ = 

(19) 

which, on use of the current conservation equation, gives the gauge invariant 
boundary condition 

A{d^F^i + ^fie^ixF^^ - Ji) = (20) 

This is seen to be consistent with the expression obtained from the field 
equations ([13|) , integrated directly over an infinitesimal region. 

3.2 Completing the action 

The gauge non-invariance of the Chern-Simons action at non-reflecting junc- 
tions is an indication of the incompleteness of the Chern-Simons theory at 
the boundary. Physically, in the region of a boundary one expects short 
wavelength modifications to come into play which will either modify the val- 
ues of physical fields or modify the action itself. To rectify the omission one 
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must either return to a more fundamental theory and rederive the correct 
effective description, or postulate the remainder of the degrees of freedom in 
such a way that gauge invariance is restored. A general requirement of the 
latter is the introduction of new variables. We now wish to discuss this latter 
problem in more detail. The general solution has been discussed by one of 
us inl 



T0[| , ignoring contact potentials. Consider the action 



S = J rf^l - If^^'F,, + ^fi{x)e,,xA'd'A^ + /(x) JM^j. (21) 



Gauge invariance implies the restriction 

1 




-{dy)e,,xd''A^ - {dj)r = (22) 



which may be written 



^{d,fi)B + l{^^^^)E' = {dJ)J' + {dJ)J\ (23) 

The conserved external source can, if necessary, be used to obtain a solv- 
able equation. If /i(x) is a function of say Xi and t then there are sufficient 
terms to find a solution for /i in terms of E"^ and B without introducing the 
source. The solution by this method will in the general case result in the 
introduction of higher derivative terms in the action. 

The role of f{x) is to act as a mediating 'leaky membrane' which separates 
the external source from the two dimensional system. Thus although is 
conserved in total, it may appear to be non-conserved via its contact with the 
two dimensional junction. Since the sole function of this term is to balance 
the gauge invariance equation, a natural boundary condition is the vanishing 
of f{x) when d^fi{x) = 0. Taking the variation to be with respect to Xi and 
t only, this is satisfied by 

/(xi,t) = ^a(9i/i) + ^/?(at/i) (24) 

for constants a and (3, hence 
(^aJ^dl - E^di + pJ^dl - Bdt)id{xi, t) + {aJ° + f3J^)didtfi{xi, t) = 0. (25) 
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This equation is too difficult to solve in the most general case, but some 
insight can be gained by noting that space and time appear symmetrically. 
Considering time-independent solutions, representing steady state solutions 
one has 

aJ\dfn) - E^din) ^ 0. (26) 
Letting y — diji and P{xi) — E'^/aJ^{xi), one has 

^ - Pi^')y - (27) 

which can be cast as a total differential by introducing the integrating factor 
exp(-/o^^ P{z)dz), 

d 



^ ,e-r-(^)- U0 (28) 



thus one has 



= ^e/;' ^(^^-^^ ^ (g{x')ei:' ^^'^A (29) 
dx Li dx y / 

where Hc and L are constants. Comparison of these equations leads to a 
differential equation for g{x^)i 

^ + ^(xW) = ^- (30) 

If, not unreasonably, E'^ and arc approximately correlated, then one may 
write P{x) = p + AP(x), for constant p and the equation for g{x^) takes the 
form 

^ + pg{x') = AP{x')g{x') (31) 
whose general solution is given by 

g{x') = r' rf;.e^(--^)(^ - AP{z)g{z)) + g^e-^^\ (32) 

When AP(a;^) = 0, this is simply the exponential decay law. In general, the 
exponential behaviour is modulated by a self-consistcntly defined function 
AP, but for any 'physical' function AP, the long term behaviour will always 
be dominated by the exponential factors. Further assumptions of symmetry 
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can be made in order to simplify even more [p!0|] . The key observation is that 
the diffusion-hke equation (26) imphes that the Chern-Simons parameter 
must fall off exponentially in space. Clearly, from the symmetry of (P3|), 
the same argument also applies for the space-independent equation for time- 
variation. The field equations in the varying region must now be solved self- 
consistently so as to record the effect of the constraint introduced through 
the solution of the gauge-invariance condition. 

Had the boundary between the regions been curved rather than linear 
then it is natural to suppose that the decay law would be modified by the 
extrinsic curvature of the interface |T6|. Since the concentration of field is 
greater around a 'corner' like feature, this would lead to an intensification of 

in this region. Although this would not seem to modify the invariance 
constraint it will, though the field equations, exacerbate the current in 
the vicinity of the interface. This increased activity must be answered by the 
sources, thus one would expect increased dissipation in such as region. 



3.3 Examples 

A simple example of the foregoing procedure can be computed in the absence 
of sources[|l^]. When /i = fi{xi, t) and E2 and B are assumed to be constants 
throughout the region of interest, the gauge invariance constraint may be 
solved together with the Bianchi identity to show that both the Chern-Simons 
coefficient and Ei are arbitrary functions of 7(i?^^xi — B~^t), for constant 7. 
This solution is extremely general and admits both longitudinal waves and 
exponential decay, but does not correspond to any obvious physical situation. 
One possibility is that it is a heuristic representation of periodic impurities 
in a quantum Hall system. 

An example of more physical relevance is the case of a time- varying Chern- 
Simons coefficient in a time-varying magnetic field. Since it is known that 
a system of spin-polarized fermions gives rise to an effective field theory 
involving a Chern-Simons term0 in the long wavelength limit, this situation 
should correspond directly to spin relaxation, or spin pumping in a 2+1 
dimensional system. Also, the construction of a model of two coupled systems 
in which spin migrates from one half-plane to the other could be composed 
of two systems of this kind, with a third junction layer of space- and time- 
varying /i describing the contact region. The constraint that / d?x^{x) be 
conserved for all times is a natural addition, implying that what leaves one 
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half-plane must end up in the other. This corresponds to the conservation of 
spin in the spin picture. Of course this is only one physical interpretation of 
such a system: the effective field equations know nothing of any microscopic 
origins and are therefore not prejudiced by the association with spin or any 
other parity violating effect. 

To solve the formal constraint, let fi = fi{B(t),t, p), p = J°, where is 



the source in (^Tf ). The constraints on the variables in the action are Bianchi 
identity 

dtB + eijdiE^ = (33) 
and the gauge invariance identity (^) 



j dV^\^[3pil-Bp^=() (34) 



where dots represent partial time derivatives and the boundary condition 
(|4D implies that 

fit) = \pp. (35) 

Eqn. (0) may be solved for p as in the preceeding subsection. It is important 
to note that the role of the gauge invariance identity is two-fold here. It is 
both an algebraic relation between variables in the action and a condition on 
the possible variations of the action. Since is an external source, one has 
that = 0, thus a variation of the action leads to 

5S = j dV,\^5A''d^F,, + ]^^5Be^'''^A^d,A^ + pe^'HA^d^Ay, 

+ \{d,p)e^''^5A^A^+^-^5B{J^A^) + f{t)J^5A}^, (36) 

where 5B = e-ijOiSAK The variation of the constraint ( |5^ yields 

The only apparent solution is 5B = 0, which implies that the vector potential 
may be written 6A^ = d^^ for some scalar function ^ and j ^ 0. The zeroth 
component A^ is unrestricted and hence there is no contradiction with (|33|). 
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Substituting this into the variation of the action p6| ) gives the field equations 
for the system. For one has, 

= JdV,i^d,E'-fi{B,t)B{t) + f{t)p^=0 (38) 

which is identical in form to the usual result for constant fi. For one has 

SS = JdV,l^{d'Oid''F,i)+fie"'\diOd.Ax 

- ^fie'^id.OAj + /(t) J^(9,0| = (39) 

Integrating by parts and using the conservation equation d^Ji = —p gives 

B + Pp = (40) 

which agrees precisely with the gauge invariance condition (|3^) up to a total 
time-derivative. Thus the consistency between the field equation and the 
gauge invariance condition is restored, by analogy with (^^. It is interest- 
ing to observe that the field equation is independent of the Chern-Simons 
coefficient, so that linearity is preserved. The equation ( p4| ) determining the 
coefficient takes the form 

fi = Cexpf^dt (41) 
J pp 

which exhibits exponential behaviour. The connection with spin relaxation 
can now be noted as follows. The variable B in the action is the effective 
electromagnetic field, not the microscopic field felt by the spins. This in- 
cludes the effect of the spin degrees of freedom. Since the coefficient of the 
Chern-Simons term is proportional in some sense to the sign of the spin and 
depends on the chemical potential of the spins p, Q this is the relevant vari- 
able to consider. One is thus interested in consistent solutions of (^Tj). As 
the effective field B tends to zero, the Chern-Simons coefficient tends to a 
constant value (typically zero). For changing B and p the external source 
is needed to drive the system. The simplest solution is for exponentially 
decaying p and B which corresponds to spin relaxation in the absence of 
an external field. If one drives the spin system with an adiabatically sinu- 
soidal time-varying magnetic field, this is refiected by an oscillatory part for 
B. This is coupled to the time- variation of p through (|38D and leads to an 
exponential lag in the response, corresponding to a hysteresis effect. 
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3.4 Contact potential: a switching junction 

A further example of physical interest arises in the case of a sharp boundary 
supporting a contact potential AAq. Taking the action {\L2^ and boundary 



condition ( ]T^ ) one has, in components 

A(5-iMo)-f = 

A{-E, + ^^^A2)-f = (42) 

and there is a step A/x at xi = 0. This system possesses certain qualities 
resembling those of a transistor or switching device: two regions of differing 
properties separated by a thin junction to which an external (bias) current 
is applied. It is straightforward to show that the picture has the properties 
of a switching device. Taking j'^ to be an external source (coupling to the 
third dimension, for instance) which acts only at the junction, conventional 
electromagnetic boundary conditions are obtained (and gauge invariance is 
restored) provided 

^AfxAo + JiAAo = -f (43) 

^Af,A, = f (44) 

where barred quantities signify the mean values of the respective parameters 
at the discontinuity {A2 = A2). Let = AAq. Since the step in the Chern- 
Simons coefficient is a physical quantity, it should be gauge invariant, thus 
equating (^) and (^) one has 

Ao Ao + doe 



and for gauge invariance 



- - ■ (46) 

A2 A2 + d29 ^ ' 

This implies a restriction on the gauge invariance of the theory at the bound- 
ary to transformations of the formQ 

e = e{A2X^ + Aot). (47) 

^For instance, 6 might be of the form exp(i(fca;^ + tot)) in which case it plays the role 
of a massless excitation. 
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This condition is satisfied by wavelike solutions, for instance. The gauge 
invariant ratio a = Aq/ A2 must be regarded as a property of a given interface 
and the gauge invariance requirement becomes 

= (48) 

The contact potential or A/i is seen to depend on the external current 
as well as the density at the boundary In particular, in the absence of 
current, the step must collapse. 

Now suppose that an external electrostatic potential V is held across this 
junction in reverse bias (opposing the potential 0). Donor charges (i.e. those 
not taken into account by the effective theory) or quasi-particles will only be 
able to surmount the potential barrier only if V" > 0. However the size of is 
controlled by the external source and thus the junction can be made to switch 
a small electric current. It is noted, on the other hand, that this 'transistor' is 
somewhat unusual since the apparent magnetic field leads also to a Hall drift 
of the electrons, and thus they drift parallel to the interface as well as across 
it. The generic behaviour described here is clearly relevant in a device which 
is populated largely by electrons of the same spin, or in an ordinary device 
with differing magnetizations in a strong magnetic field. It is known that, in 
a system of fermions with predominantly one spin direction, a Chern-Simons 
term is induced by radiative corrections in the long wavelength limit, the 
coefficient of which is proportional to the two-dimensional spin eigenvalue. 
The above scenario then describes a spin-switch, i.e. a semi-conductor-like 
switch which is controlled entirely by parity-violating spin effects. Some 
experimental evidence for this exists already. 

The flow in and out of the system by sources at the boundary, in our 
notation, has the appearance of a charged current. This is because is 
formally conserved. However, the variable coupling f{x) implies that is 
not conserved on the boundary. There, the lack of manifest gauge invariance 
implies that the compensating current could appear neutral, as seen from 
the junction's perspective. Thus need not be interpreted as a current 
of electrons. The dissipation could relate to radiated energy generated by 
spin flip transitions in moving from one region to the other. If an electric 
current is prevented from flowing freely then it would be expected, from the 
preceding discussion, that there would be some resistance to the transport 
across the junction, since a tendency toward a reflective boundary condition 
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would arise. This could lead to a component of Ji arising, implying through 
(p5D an accumulation of spins at the junction. While these would decay by 
diffusion eventually, there is the possibility of some hysteresis depending on 
the conductivity of the material at the boundary. 

Finally, it is important to bear in mind that the present theoretical model 
exists only at the level of an effective field theory, not a microscopic one. The 
motion of dressed quasi-particles through an interface of changing physical 
properties must take place through the interaction with some third party, 
since the particles must be 'undressed' and 'redressed' in order to make the 
transition across the barrier. In the above construction, this occurs through 
an interaction with the source. The source itself consists formally of quasi- 
particles characterized by the mediating value of the physical parameter Ji. 
In the above discussion we have adopted the classical viewpoint of an effective 
field theory to avoid dealing with the subtleties of quantum tunneling, but 
in the final analysis one is interested in dispensing with the effective theory 
and finding the appropriate microscopic one. Tunneling processes will be 
relevant at this level. 



3.5 Experimental evidence 

After predicting the diffusive junction behaviour and the above switch model, 
we learned of two separate experiments which have a direct relevance to these 
findings. Experiments by Kane et al.[^] consider the interface between two 
quantum Hall systems with different filling fractions^. This corresponds to 
differing Chern-Simons coefficients according to the arguments in refs. [^, ^]. 
The key observation is that the aligned spins break parity invariance and 
that the magnitude of this breaking depends on the number of spins[0. Here 
one observes a diode-like behaviour, as expected from the general arguments 
above. It is interesting that this experiment shows evidence for a polarization 
of nuclear spins at the junction. This coupling to nuclear spin acts as source 
(sink) of the polarization in crossing the junction and constitutes the 'external 
system' in the preceding source language]^, |IU . 



A second set of experiments by Johnson relates more directly to the 
switching phenomenon!]^. Here a device is constructed in which spin-polarized 
charge carriers inhabit a thin gold layer with a two dimensional symmetry. 



^We thank B. Halpcrin and T. Finstad for pointing out these references. 
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Spins subsequently migrate across a junction into a region of different or op- 
posite spin||] by passing through a finite width junction which is coupled to 
an external source of current. A 'spin bottleneck' phenomenon is observed 
which prevents the non-equilibrium junction from decaying when the source 
is active and switching is indeed observed, in accordance with the above 
discussion. 



4 Non-minimal coupling 

In this section we study a related system which avoids some of the gauge 
invariance problems of the Chern-Simons system. 

4.1 Gauge Invariance 

Consider the action 

S = J rf^l - If^^F^, + ^K{x)e>^'''j^F,^^ (49) 

for some conserved current J'^. This model has been studied in a num- 
ber of recent works [p!7| , |l^, |l^, ^ The coupling of the current to the 
dual of the field strength leads to a local magnetic moment type interaction 
for scalar particles[0]. In 2-1-1 dimensions this leads to an induced phase 



analogous to the anyon phase, known as the Aharonov-Casher phase [PH 
An interesting duality exists between this term and the usual Chern-Simons 
type parity breaking term. If one begins with scalar electrodynamics coupled 
to a Chern-Simons term, the non-minimal coupling is induced by radiative 



corrections [|T8[. The coefficient is related to the inverse of the Chern-Simons 
coefficient and the electric charge and does not vanish in the long wavelength 
limit. Similarly, if one begins with the non- minimal coupling and computes 
one-loop corrections, a regular Chern-Simons term is induced. Apart from 
these properties, this non-minimal coupling is interesting in the present con- 
text due to its manifest gauge invariance, even in the presence of a boundary. 
It is natural to compare the boundary properties of the above model with 
the more usual Chern-Simons term. 

Applying the action principle to (HOf), one has 



d^lF^"^ + K{x)e>"'^J^] = (50) 
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and associated boundary condition, for a junction at xi = 0, 



A(F^-' + K{x)e"'^J,) = 0. (51) 



It is now possible to compare the action ( ^Uj ) to the Chern-Simons action 
coupled to a source (0). A direct comparison is only possible if the current 
itself includes terms involving the gauge field. This is the case for a complex 
scalar field, for example. Then one has 

J^ = z[$(D^$)t-$t(D^$)], (52) 

where Z?^ = — ieA^, or in the unitary parameterization $ = pe*^, 

J, = 2p\d,e)~eA^]. (53) 

The non-minimal term in ( ^9]) then takes the form 

2K{x)p'F*^'{d^e - eA^) (54) 

where F*^^ = ^e^^'^F^x. Comparing this to (^) and derived quantities, it is 
possible to identify 

/i(x) = -AeK{x)p'^ (55) 
/(x)J^ = -2e^^\d,t,p){d,9) (56) 

which indicates that the non-minimal coupling can be regarded as a normal 
Chern-Simons term with a variable coefficient plus an external massless scalar 
field 6 which lives in the region V/i(a;) 7^ (for example on boundaries). This 
is consistent with the discussion in section 3.2 and in refs. P, ITOl. 



4.2 Atomic spring model 

For the purpose of determining bound state properties for comparison with 
ref . [|TT| , it is interesting to consider a toy model of a two dimensional medium 
in which electrons are bound to their parent atoms by means of a pseudo- 
harmonic potential. If s* is the displacement vector of a single electron then 
under a Lorentz force, 

{dl + ^dt + ul)s' = --{E' + e,,BAs') (57) 
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where is a constant external magnetic field, uOq = k/m for spring constant 
k and 7 is a damping factor. The current arising from this motion can be 
characterized by 

where p is the total charge density (which is usually zero, including the effect 
of background charges) and is the number of optically active electrons. 
Adding a minimal gauge coupling J^A^ to (pQ]), and omitting dielectric po- 
larization effects one has the field equation for constant k 

d^F^"^ + Ke>"'^d^J^ = (59) 



and boundary condition given by (pll). We shall use this model below in 
order to determine the optical properties of this system. Such properties 



were previously calculated for the usual Chern-Simons term in ref. |11 



4.3 Complex Scalar Field 

A more realistic quantizable field theory can be obtained by considering a 
scalar field coupled minimally to a gauge field. This is closely related to the 
effective P and T breaking theory introduced in ref. ||22| for superconduc- 
tivity, and the to the Landau- Ginsburg theory of the Hall effect [^. The 
scalar field represents the collective excitations of an unknown microscopic 
system in the long wavelength limit. In the Chern-Simons case, this model 
has already been shown to lead to a Faraday rotation in reflected waves p^] . 
The principal difference from the atomic spring model is that the scalar field 
is a superconductor-insulator. There is no in-built mechanism for dissipation 
at zero temperature, hence no finite conductivity. Moreover, the model does 
not describe bound states. In the superconductive regime, waves will not 
propagate inside the material, but extra-planar waves can be reflected off 
a two-dimensional surface in which the scalar field lives. In the insulating 
regime, wavelike-solutions can only exist under special circumstances which 
will be described below. 

The action for a complex scalar field is defined by 

S = J dV,i^{D^^<^y{D^^) - m2$t$ _ ^($t$)2 _ If'^-f^, + ^Ke^-^J^F,^^ 

(60) 
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where the current is given by (0) and = <9^ — ieA^. It is noteworthy 
that the matter couphng involves a quadratic dependence on the gauge field 
that cannot be written strictly in the form J^A^. The field equation for the 
scalar field is given by 

(D^ + m^)^ + 2iKF*^{D^,^) + -(^t^)"!) + = (61) 

3 

with associated boundary condition at Xi = 0, 

A(Di + l«:ei,;,F^^)$ = 0. (62) 

For the gauge field, one obtains 

d^F^"^ + Ke'^'^df.J^ - 2eKF*^$^$ = eJ^ (63) 

provided 

A{F^^ + Ke^'^J,) = 0. (64) 



5 Energy Momentum Tensor 

The gauge non-invariance of the Chern-Simons Lagrangian at a boundary 
is accompanied by a discontinuity in the electromagnetic Poynting vector 
AS*! = A(i?2-B)|^ H^- This is another indication that the theory is not 
complete at the boundary. To properly understand the energy fiow in the 
present models one should supplement the usual electromagnetic Poynting 
flow by a contribution from the parity violating term. A consideration of the 
energy-momentum tensor as defined through Nother's theorem leads to the 
relevant conserved quantitie^. The generalized force law for the system is 
given by 

r = J dv^dj-"^' (65) 

which, in the case of the electromagnetic field, gives rise to the Lorentz force 
law = — / dVxJ^F^^, where dv^ is a spatial volume element. Any anoma- 
lies in 9^^ could give rise to corrections to this force law and are therefore 

■^If one defines the energy-momentum tensor by T^" — for the metric g^^, then 

there is no contribution from parity violating terms, since e^"'*' transforms hke a tensor 
density and is subsequently independent of the metric. However, this should be regarded 
as a failure of the variational definition rather than a reason to disregard the extra terms. 
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important to the discussion of boundary effects. A generic feature of the 
energy-momentum tensor in the present models is that it is non-symmetric 
and gauge-dependent. This signals a breakdown of Lorentz invariance as well 
as gauge invariance. That these two principles should be broken simultane- 
ously is reminiscent of the deficiencies of the canonical energy momentum 
tensor |]12|, |l^ and thus in what follows we shall adopt the covariant proce- 



dure of refs. |T2| and |]T3| to obtain directly the normally symmetric 'Belifante 
energy-momentum tensor'. 

To construct 9^^ , one notes that the covariant vector potential trans- 
forms like a 3-vector only up to a gauge transformation. Nother's theorem 
implies that under a Lorentz transformation 

_^ ^f^' = + Sx^" (66) 
A^'ix) A^"{x) = A^'ix) + 6A^{x) (67) 

a symmetric theory is characterized by the continuity equation d^C^ = 0, 
where 

where S = J dVxC The defining equation for 6^'^ is then 

= -e'^'Sx^. (69) 
The gauge-Lorentz invariant restriction 6A^ = F^"5xa implies that 

If a boundary is introduced then Lorentz invariance is explicitly violated and 
this expression loses its immediate interpretation as a consequence of the lack 
of translational invariance. In the absence of a parity violating term one has 
the usual electromagnetic energy momentum tensor 

e^^,^ = F-^F^^ + -^r^^^FP^F^,. (71) 

The corrections A^^*^ to this due to the parity violating terms will now be 
discussed. 
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5.1 Chern- Simons 



Use of the Chern-Simons Lagrangian in ([70| ) gives 

A^A- = - ^fxri'^-'eP'^'^ApF,,. (72) 

This contribution is not symmetric, but, using the Bianchi identity d(^^FyX) = 
e'^'^^d^Fyx = it is seen to be gauge-invariant provided the coefficient /i is 
a constant. The integral of the divergence of the full 9^^ = 9^^^^ + A6'^'^ is 
conserved in the absence of sources, 



/ 



dV^d^9^'^ = (73) 



where the field equations (|T^) have been used. Adding sources produces the 
conventional modification 

j dV^d^9^^ = - j dVJ^F-^ (74) 

which is the standard Lorentz force law. The latter result implies that no 
modification of the force law is required in Maxwell-Chern-Simons theory, 
except at boundaries or in regions of variable /i(x) in which case one finds 
that 

j dV^d,9^'' = j rfK|^(x)F-(^n + \{d^A^))^. (75) 

Under a gauge transformation ^ + d^Xi this term changes by the 
gradient of — |((9°"/i(x))F^X5 which may be compared to (p2|) for zero source. 
This identification shows explicitly that the lack of energy conservation at 
the boundary is related to the lack of gauge invariance, as previously argued 
in ref. and leads to a generalized force on the system. 



In ref. the presence of an instability in the dispersion was found. 
Although some evidence for this can be seen in 9^^ (which has no definite 
sign), this does not cause any problems in the present work. 

5.2 Non-minimal term 

For the non-minimal Lagrangian, one has 

A^/^- = ne'^-'^J^F"^ - ^r/^^/te^'^V^F^A, (76) 
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and it is assumed that is independent of A^^^. This correction is non- 
symmetric, but is exphcitly gauge invariant. The integral over the divergence 
of the total energy-momentum tensor is conserved, provided k is constant, 
and the addition of sources leads to the usual Lorentz modification (ff^). For 
non-constant one has 

j dV^d^e^'' = -J dK.«:e^'^'^|(5^Jp)n - lid''Jp)F,,'^. (77) 

This term is gauge invariant and represents a modification of the standard 
Lorentz force wherever k is a function of x^. 

In ref. it was shown that a Casimir force is necessarily present on a 
gauge-invariant reflective boundary due to quantum or thermal fluctuations 
in Maxwell-Chern-Simons theory. The force arises because there is a mass 
or gap-mismatch in the spectra for the two sides of the junction. Here the 
spectrum of the non-minimal coupling is notably gapless (massless) in the 
absence of an ordered phase (spontaneous symmetry breaking). See equation 



(|106|) . This seems also to be confirmed by refs. |T9|, ^ at one loop, but 
here one notes that scalar and gauge loops are inextricably linked by the 
non-minimal term so that it is not possible to consider the gauge field in 
isolation. Further investigations would be required to determine the absence 
of a gap due to fluctuations in general. 

Finally, 6*°° does not have a definite sign, which suggests the possibility 
of an instability for certain values of the sources. No instability is found in 



the dispersion ( 106 , 120 ) however. In the case of normally impinging waves 



for inhomogeneous currents, this is less clear 



6 Electromagnetic waves 



In an earlier paper |TT|] one of us has considered the properties of wavelike 
solutions of the Lagrangian ([T^). Related work has also been carried out in 
3-1-1 dimensions ||14||. There it was remarked that a future problem would 
be to consider the analogue of Fresnel's equations at a material interface. 
Since the naive boundary condition for the Chern-Simons model is not gauge 
invariant, it is not immediately clear how to proceed. To satisfy the gauge 
invariant boundary condition (p!7[) , there must be total reflection from the 
line Xi = 0. If on the other hand one couples to an external source as in (|21|), 
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then the gauge non-invariant parts of the boundary condition vanish and one 
is left with the normal electromagnetic boundary conditions. The virtue of 
the non-minimal model is its automatic gauge invariance at the boundary. 



6.1 Maxwell- C hern- Simons atomic spring model 

We begin with the unmodified Chern-Simons theory and consider the case 
in which the coefficient is constant. 

6.1.1 Plancir dispersion relation 

The dispersion relation is obtained on substituting the trial solution 

E' = £;^e*(^'^'-'^*) (78) 
B = 5oe*('=^^^-'"*) + 5c (79) 

Ji ^ ji^i{kiXi-U,t) 

P = Poe^^'^'^^-'^^ + Pc (81) 
into the components of the field equations 

Be = --Pc (82) 
P 

ikEow = pBo + po (83) 

iuEow - nEo_i_ = Jo|| (84) 

-ikBo + iujEo± + /iEqii = Jo± (85) 

where the parallel and perpendicular projections are defined through the 
relations 

k'Ei, = kEow 
eijk'El = kE^^ (86) 

and similar ones for the current. If one neglects the oscillatory part of the 
magnetic field from the equations of motion (which implies that the magnet- 
ically induced current-response is small) then these equations, together with 
(57) and (58), can be manipulated so as to eliminate all variables except the 
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electric field, which then satisfies the equation 




/ iuomW + iNe^uj — i^uoeBc —fimW — eBJi^ + uPeB^ 

\ mW - uj'^eBc + iuj)mW + iNe^uj - ifiooeBc 

where W = (— a;^ — iju + uJq). Demanding the vanishing of the determinant 
of this matrix leads immediately to the dispersion relation 

(88) 

where Up = Ne'^/m and ujl = eBc/m. This result has been given in ref. 
[H]. The refractive index is given hj n = k/uj. The real and imaginary 



parts of the refractive index are plotted against uo for various values of in 
Figs, la and lb. At a junction between two regions, this dispersion relation 
applies provided the interface is sharp. Sources are then needed to balance 
the requirements of gauge invariance, but their introduction leads only to 
normal electromagnetic boundary conditions, thus there is no modification 
to the laws of refiection or refraction and Fresnel's relations are given by the 
usual formulas (see below). Thus when two dimensional waves strike the 
interface, currents are set up along the interface. These currents must either 
disperse into the two dimensional system or out into the external system. If 
no sources are introduced, the dispersion is only given by (|8^) at t = 0. The 
subsequent decay of the interface then modifies the dispersion in a non-linear 
way until a situation of equilibrium is reached. 



6.1.2 Extra-planar dispersion relation 

To show that the Chern-Simons term gives rise to a modified Faraday effectQ, 
one embeds the planar model into a three dimensional space and directs plane 
polarized waves so that they impinge normally to the plane. Since the matter 
fields describe an ostensibly two dimensional system all currents are restricted 
to the plane. The embedded action is given by 

S = j dV,\^ - \f^'F^, - JM, + ^/ie^,),3^^^''} (89) 

^This is also referred to as optical activity or circular birefringence. 
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where careted indices run over the third spatial dimension in addition to their 
usual values and p, has dimension one. The components of the field equation 
are then given by 

d-E^-jlB' = p 

d,B^ + e^fPiE^ = 

d-B^ = 0. (90) 

Next consider a solution of the form 

E' = Ele''^^'-'''^ (91) 

with corresponding expressions for the magnetic field and current. To obtain 
the dispersion relation between k and uj it is sufficient to consider the planar 
components of the electric field. It is convenient to define complex variables 
E = E^+jE"^ where = —1, but ij = ji ^ — iQand corresponding variables 



for the displacement vector s|15]. With these variables the equation of motion 



for the electrons ( p7D becomes (neglecting damping terms) 

{-u^ + ul)s = -^{E- ]B,dts). (92) 

Combining the second and third equations in (pO|) gives 

Neu^s = {uj^ - k^)E - ijujjjE. (93) 

Finally, eliminating one of the variables from the last two equations gives the 
dispersion relation 

k =u; + — — ^-2j(- + ^„ 



^'(-Cu2 + ^2)2_^2^2 -^V^ ^'(-C^2 + ^2)2_^2^2^j- 

(94) 

The real and imaginary parts of the refractive index n = k/uj are plotted 
against uj for various values of jl in Figs. 2a and 2b. The combination ij 

^Note that if these imaginary units are treated as anticommuting numbers they re- 
produce the su(2) algebra of rotations. Here, as commuting numbers they represent a 
reahzation of the group U{1) x U{\). 
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guarantees energy conservation and implies that the wavenumber has the 
form 

which shows that the transmitted wave suffers a rotation of its plane of 
polarization 

E = Eoe'^^'-'-'^'\cos kijz + j sin kijz). (96) 

It is normal to define the measurable rotation angle in terms of 'Verdet's 
constant' V by 

Xrot = VB,L (97) 

where L is the distance travelled through a stack of planar systems, but since 
the Chern-Simons term results in such a rotation even in the absence of Be, 
it is sufficient to consider the angle rotated per unit length kij. 



6.2 Non-minimal atomic spring model 
6.2.1 Planar dispersion relation 

Wavelike solutions to the non-minimal action (^9|) are obtained by the same 



method as in ref. [11| . It is curious to note that the presence of the curl of the 
current makes the field equations sensitive to anisotropy and inhomogeneities 
in the the matter field. The components of (p9|) give, 



diE' - Keijdi.P = p (98) 
-dtE' + eijdjB + KEi, {djp + dtJ') = f (99) 

Using the wave ansatz ([r^)-(|81|), and considering only the oscillating parts 
leads to 

ikE^w — iKkJo± = po (100) 
iuEQW — iu!KJoj_ = Jo|| (101) 
iujEoj_ — ikBo — inkpo + iKuJo\\ = Jo± (102) 
ikEo± = iuBo (103) 

where the last equation follows from the (|33|) . The continuity equation for 
the current (which is not an independent equation, but follows from ( |100| ) 
and ( |101| ) above) may be written 

kJp = iupo. (104) 
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Eliminating the electric field in favour of the current, one obtains after some 
effort the matrix equation 

/ mW + iVe^ iNe^UK - iuoeB^ \ ( 'hw \ _ 
\ lueB.-iNe^cuK mW + Ne^J^ ) \ J^l ) ^ 

The vanishing of the determinant of this matrix yields the dispersion relation 



k'=UJ 



1 + 



(106) 



The real and imaginary parts of the refractive index n — k/uj are plotted 
against u for various values of n in Figs. 3a and 3b. 
It is noteworthy that, on making the definitions 

H = B + Kp (107) 

the wave equation for the transverse magnetic waves takes the usual form 

{d^ - V^)H = SijdiJ', (108) 

and the Bianchi identity is unchanged, implying that the normal laws of 
reflection and refraction apply: 

sin^t rii 

^ = - 109 

sm Oi n2 

9i = 9r (110) 

where the angles 6 refer to the incident, transmitted and reflected angles 
to the normal. It is quite possible that a non-linear scalar field theory (for 
instance, with a term or higher power) would violate this law, if wave 
solutions can even by found. 

Prom the linearity, it is possible to write 

{Eoi - Eor) COS 9i = Eot cos 9t (111) 

Thus, defining the wave impedance Z = E / H = Hr/n^, where /irH = B, 
and noting that the boundary condition may be written AH — 0, one has 
the familiar result that 

^ + ^-^. (112) 

Z/i Zj\ Zj2 
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Combining ( |111| ) and ( |112| ) gives Fresnels standard relations for the trans- 
mission and reflection coefficients 



Eor 

Eoi 
Eoi 

6.2.2 Extra-planar dispersion relation 

To compute the Faraday effect, the planar theory is embedded in a three 
dimensional space, with waves impinging normally, as before. The action is 
therefore 

S = l rf^l - iF^-'Ef^o - J'' A, + he^.Xs-J'F'''^ (US) 

where careted indices include the third dimension. The currents are assumed 
to lie purely in the plane. The components of the field equation and Bianchi 
identity are now given by 

diE"- — keijdiP = p 
-dtE^ + eiy^d-B^ + ke^jidjp + dtf) = f 

dtB^ + e^idiE^ = 

d-B' = 0. (116) 

The oscillating magnetic field can be eliminated to yield a wave equation for 
the electric field 

^ -^i^^iikz-ujt) ^ (117) 

given by 

i-d^ + V^)E' + keijidtdjp + dp') = dtf. (118) 

It is pertinent to note that the presence of spatial derivatives of the current 
in the above expression implies that the system is sensitive to inhomgeneities 
and anisotropy. A proper description of anisotropy cannot be obtained from 
the foregoing equations, since the spring constant k = oj^m would need to be 
different in the xi and X2 directions to make the assumption self-consistent. 



ZicosOi — Z2 cos Ot 
Z2 cos 6t + Zi cos 6i 

2Z2 cosOi 
Z2 cos 9t + Zi cos 9i 



(113) 

(114) 
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6.2.3 Homogeneous isotropic medium 

In a homogenous, isotropic medium, the spatial derivatives of the current 
vanish identically. Using the complex coordinate representation introduced 
earlier one has 

(tu^ - k'^)E + ijRuj^Nes = u^Nes. (119) 



Using the equation of motion for s (92), it is possible to eliminate E giving 
immediately the dispersion relation 



(-a;2 + ijuLUJ + uJq] 



(120) 



or 



^' = -'{^^^ \\ (121) 



The real and imaginary parts of the refractive index n = k/uj are plotted 
against uj for various values of k in Figs. 4a and 4b. 

kij is the rotation per unit length through a layered system, where k = 

6.2.4 Homogeneous, anisotropic medium 

If the wavelength of waves is small compared to the inhomogeneities of the 
material medium, the dispersion becomes sensitized to the structure. Con- 
sider the case of a uni-axial crystal which permits inhomogenities of the cur- 
rent in a preferred direction Xi. Using the trial solution E = i?oexp(z(fc^z + 
kiXi — ut)) one easily obtains 



2 kl - kl)E^ - iNcKuj^s^ = Neuj^s^ 



kl-kl)E^ + iukNe{uj^ -kl)s^ = Neuh^. (122) 



The complex method is not appropriate here, owing to the lack of symmetry. 
Nevertheless, it is possible to eliminate the electric field component-wise and 
solve the determinental equation for the matrix coefficient of and s^, giving 
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where a = uJq — uj"^ — cu^culk, b = — 2a;pa;^((a;Q — uj'^) + kujii^i^ — uj)) and 
c = ujpUJ^il + k^ikl — oj"^))- The dispersion continues to exhibit birefringence, 
but now with a marked asymmetry. The ki terms either reduce or enhance 
the electron mobihty, depending on the sign of the magnetic field uj^. In 
particular, it is seen that ki acts as an effective mass gap for the dispersion. 



6.2.5 Ohmic conductors 

The conducting limit of the previous results could in principle be obtained 
from the uOq^{) limit of the atomic spring model. It is useful to reexpress 
the result in terms of the more familiar conductivity. Let the anisotropic 
conductivity tensor be defined by its projected components 



J\\ 



ail Ell. 



(124) 



then from (|100D-(|103D one has 




-Ka± 

iuja 




(125) 



giving the dispersion relation 



1 + 



-K 



(126) 



which reduces to standard results on setting n = 0, a\\ = 0. What is inter- 
esting here is that the longitudinal current plays a role in the dispersion. In 
the vicinity of a boundary, like the edge of a finite sample, the simple split 
into (T|| and a± must break down. Close to the edge, the longitudinal con- 
ductivity must tend toward zero and be replaced by an enhanced transverse 
conductivity. This corresponds to a modulation of the charge at the edge of 
the sample, which then spreads out to form surface density waves. 

The Faraday effect is straightforwardly obtained from ( 119 ) on substi- 
tuting J = aE] we ignore the role of anisotropy and inhomogeneities here. 
Then, straightforwardly 



UJ 



1 + jKa + i- 



(127) 
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which indicates that the rotation of the polarization plane is accompanied 
by damping and reflection. It is interesting to note that the combination jn 
implies that the reflective properties are unaffected by k = whereas the 
polarization effect is entirely due to k 7^ 0. 



6.3 Non-minimal Complex scalar field 

Wavelike solutions for the electromagnetic field are not a general feature of 
the non-minimally coupled scalar field theory. Consider the case in which the 
collective field mode is covariantly constant i.e. -D^$ = 0. Then the scalar 
field equation implies that 

$2 = f; + m^). (128) 

A 

If the field strength is oscillating, this makes most sense when d^n = 0. In a 
superconducting phase there is clearly no wave propagation in 2 dimensions. 
There is a regime however in which propagating solutions can be obtained. 
For covariantly constant $ the current vanishes but is an invariant con- 
stant. The field equations for the electromagnetic field are then 

d^F^^ - 2eKF*^$T<l> = 0. (129) 

Taking the spatial components of this equation together with the Bianchi 
identity ( P5D leads to a determinental equation for the dispersion relation 



which is notably similar to that for the superconductor model in ref. 
indeed it is regular Chern-Simons dispersion for /x = 2eK$'^$. 



7 Discussion 

We have examined some of the consequences of parity violation in the vicinity 
of junctions and boundaries. Using the principles of gauge invariance and 
energy conservation, we derive the acceptable behaviour of two models for an 
effective P and T breaking theory. In the case of regular the Chern-Simons 
term, it is found that dissipation must be a feature close to a boundary. This 
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dissipation can either be a destructive dissipation - that is, one which tends 
to erode the boundary itself, or a stable transfer of current to an external 
system by mediating sources. It is interesting to note that, in the quantum 
Hall system, edge currents can be measured and that dissipation is observed 
at the boundaries of the sample, concentrated at the corners, which is at least 
in qualitative agreement with the picture conveyed here. The non-minimal 
coupling has by nature edge currents, but these are related unambiguously 
to the magnetic moment coupling and result in no dissipation. On the other 
hand, the presence of a 'curl' of the current implies a mixing of charge at the 
boundary with the bulk charge. 

The diffusion-like behaviour of a regular Chern-Simons interface, together 
with the need for an external source of current suggests a transistor like 
behaviour. While this behaviour is in itself amusing, since it is derived from 
a consistency argument in an effective theory, it appears to have relevance to 
experimental devices in which the charge carriers in different isolated regions 
have predominantly the same spin^, |^. An electron passing from one region 
to another must then flip spin, requiring either an input or a drain of energy. 

We have considered the effect of penetration of material samples by elec- 
tromagentic radiation. Waves inside P and T breaking media are no longer 
transverse, but no essential modification of the usual laws of reflection or 
refraction is noted. Dispersion is qualitatively different for the two models 
considered. In the regular Chern-Simons theory, the Chern-Simons param- 
eter appears mainly as a mass or gap term in the dispersion relation. The 
relation for the non-minimal model does not appear to possess such a gap. 
Both models exhibit the required optical activity, or circular birefringence. 

Finally we note that similar studies of boundary effects in parity violat- 
ing models have been made in refs |2^, and ||2^. In the former case it is 
shown that a connection exists between the chiral anomaly in 2n -f 2 dimen- 
sions and the Chern-Simons gauge anomaly in 2n + 1 dimensions, at least 
to first order in a derivative expansion; when chiral models are calculated 
non-pertubatively, one also finds that the Chern-Simons form is modified by 
higher order terms and becomes the 77-invariant. In the latter case, edge 
currents are argued by appealing to linear response theory and gauge invari- 
ance. In both of these cases, currents are responsible for balancing the gauge 
invariance constraints. In this paper further solutions are found which do not 
require the inclusion of additional currents and a new physical interpretation 
is given to the gauge invariance problem. 
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Figure Captions 

la. The real part of the refractive index n = k/uj as determined from Eqn. 
(88) plotted against uj for various values of /x. We use luq = 1, cUp = .Iluq, 
cul — -OIcuq and 7 = .01a;o. 

lb. The imaginary part of the refractive index n = k/iv as determined from 
Eqn. (88) plotted against cu for various values of /i. We use ljq = 1, 
cUp — .Icuo, ujl — .QIujq and 7 = .01a;o. 
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2a. The real part of the refractive index n = k/oo as determined from Eqn. 
(94) plotted against oj for various values of ji. We use uiq = 1, cUp = .lo^o, 
ujl — .QIujq and 7 = 0. 

2b. The imaginary part of the refractive index n — k/uj as determined from 
Eqn. (94) plotted against uj for various values of pL. We use ujq — 1, 
ujp — .Iujq, ujl — .01a;o and 7 = 0. 

3a. The real part of the refractive index n — k/uj as determined from Eqn. 
(106) plotted against u for various values of k. We use ujq — 1, 
ujp — .lujQ, ujl — .01a;o and 7 = .QIujq. 

3b. The imaginary part of the refractive index n = k/oo as determined from 
Eqn. (106) plotted against uj for various values of k. We use Uq — 1, 
ujp = .Iujq, ujl = -01(^0 and 7 = .01a;o- 

4a. The real part of the refractive index n = k/oo as determined from Eqn. 
(121) plotted against uo for various values of k. We use uoq = 1, 
oOp — .looo, ool — .01a;o and 7 = 0. 

4b. The imaginary part of the refractive index n = k/oo as determined from 
Eqn. (121) plotted against 00 for various values of k. We use ooq — 1, 
OOp — .looo, ool — .01a;o and 7 = .0. 
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